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Abstract. This paper describes the algorithm and examines the performance of an 
IMRT beam-angle optimization (BAO) system. In this algorithm successive sets of 
beam angles are selected from a set of predefined directions using a fast simulated 
annealing (FSA) algorithm. An IMRT beam-profile optimization is performed on each 
generated set of beams. The IMRT optimization is accelerated by using a fast dose 
calculation method that utilizes a precomputed dose kernel. A compact kernel is 
constructed for each of the predefined beams prior to starting the FSA algorithm. 
The IMRT optimizations during the BAO are then performed using these kernels in 
a fast dose calculation engine. This technique allows the IMRT optimization to be 
performed more than two orders of magnitude faster than a similar optimization that 
uses a convolution dose calculation engine. Any type of optimization criteria present 
in the IMRT system can be used in this BAO system. An objective function based on 
clinically-relevant dose-volume (DV) criteria is used in this study. This facilitates the 
comparison between a BAO plan and the corresponding plan produced by a planner 
since the latter is usually optimized using a DV-based objective function. A simple 
prostate case and a complex head-and-neck (HN) case were used to evaluate the 
usefulness and performance of this BAO method. For the prostate case we compared 
the BAO results for 3, 5, and 7 coplanar beams with that of the same number of 
equispaced coplanar beams. For the HN case we compare the BAO results for 7 
and 9 non-coplanar beams with that of 9 equispaced coplanar beams. In each case 
the BAO algorithm was allowed to search up to 1000 different sets of beams. The 
BAO for the prostate cases were finished in about 1-2 hours on a moderate 400 MHz 
workstation while that for the head-and-neck cases were completed in 13-17 hours on 
a 750 MHz machine. No a priori beam-selection criteria have been used in achieving 
this performance. In both the prostate and the head-and-neck cases, BAO is shown to 
provide improvements in plan quality over that of the equispaced beams. The use of 
DV-based objective function also allows us to study the dependence of the improvement 
of plan quality offered by BAO on the DV criteria used in the optimization. We found 
that BAO is especially useful for cases that require strong DV criteria. The main 
advantages of this BAO system are its speed and its direct link to a clinical IMRT 
system. 
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Beam Angle Optimization 


CRT 


Conformal Radiotherapy 


CSA 


Classical Simulated Annealing 


DV 


Dose Volume 


DVH 


Dose Volume Histogram 
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Equivalent Uniform Dose 
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l^ast Simulated Annealing 


HN 


Head and Neck 


IMRT 


Intensity Modulated Radiation Therapy 


lO-SC 


intensity Optimization with Superposition Convolution 


lO-TLP 


Intensity Optimization with Table Lookup 


OAR 


Organ at Risk 


OF 


Objective Function 


SA 


Simulated Annealing 


SC 


Superposition Convolution 


TLP 


Table Lookup 


TPS 


Treatment Planning System 


VOI 


Volume of Interest 



1. Introduction 

The beam-angle optimization (BAO) in intensity-modulated radiation therapy (IMRT) 
aims to find, for each patient-specific case, a set of beam angles that can produce a three- 
dimensional (3D) dose distribution that conforms most closely to the prescribed one. 
In general, selecting the best set of beam angles that can achieve this goal is a highly 
non- intuitive process. This is already true in 3D conformal radiation therapy (CRT), 
where each beam has a preset profile inside its beam's-eye-view aperture and only the 
beam weights are optimized, but especially more so in IMRT. The optimal distribution 
of the total dose in IMRT is the result of a complex interplay among a large number 
of beamlet intensities from multiple beam angles and therefore there is no a priori way 
of predicting the dependence of the objective function (OF) on the beam angles. The 
optimizer has to balance the requirement of creating a prescribed dose distribution at 
the target with that of giving as little dose as possible to the surrounding normal tissue 
and organs at risk (OARs). Stein et al (1997) provided a well-known counter-intuitive 
example of the necessary compromise that often has to be struck to achieve this balance 
between target coverage and OAR sparing: sometimes the beams have to go directly 
through an OAR to achieve an optimal distribution. 

The complex dependence of the dose distribution on the beam angles and their 
beamlet intensities gives rise to a multiple local minima problem. The fact that the 
OF has a highly non-convex dependence on its variables has so far precluded any 
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deterministic approach to its minimization and necessitated a brute-force assail to this 
problem, namely a global search method. This approach requires one to perform an 
IMRT optimization on each one of a large number of possible combinations of beam 
angles. Two main hurdles have to be overcome, however, in making this highly-iterative 
approach practical and can be used routinely in a clinical setting. The first hurdle 
is the substantial amount of time that is required to perform an IMRT beam-profile 
optimization for a single set of beam angles. Depending on the complexity of the case 
and the algorithm used in the optimization, the required time already ranges from 
several minutes to a few hours on a typical workstation. The second hurdle is the huge 
size of search space involved in the beam-angle search. To rigorously search for, e.g., 
the best 5 gantry angles out of 180 possible directions would require one to search over 
Cfgo(~ 10^) possible combinations. A BAO system that is based on an exhaustive search 
would therefore, most likely, be too slow for utilization in a clinical setting. 

The recent past has seen a limited but steady flow of publications on this topic. 
A few papers have provided theoretical analysis of the BAO problem (Bortfeld and 
Schlegel 1993, Jackson et al 1994, Braunstein and Levine 2000). Bortfeld and Schlegel 
(1993) considered the BAO for the 3D CRT, as opposed to IMRT, and concluded that, 
for dose-based criteria, the optimum beam configuration for multiple-beam irradiation 
tends to be an even distribution over the 27r angular range for the gantry. 

Rowbottom et al (1998, 1999a, 1999b) have proposed several approaches for BAO as 
applied to the 3D CRT. Dose-based OF was used in the optimization of the beam weights 
and simulated annealing algorithm was used to select the beam orientations in a reduced 
search space of coplanar (Rowbottom et al 1998) and non-coplanar (Rowbottom et al 
1999a) beams. Some a priori criteria were also used in the beam-orientation selection: 
the beam orientations were constrained to be separated by at least 35° and parallel- 
opposed beams were disallowed. In addition, a single-beam cost function (Oldham 1998), 
based largely on the geometry of the target, was used to rank the beam orientations in 
the search space according to their likelihood for being a member of the optimal set and 
only the top 10% of the set were searched in the BAO. Rowbottom et al (1999b) also 
attempted a more intelligent search using an artificial neural network technique. The 
results of this search, however, seem to be inferior to those that were obtained using 
their earlier simulated annealing method. 

Pugachev and Xing recently proposed several approaches to BAO for the IMRT 
(Pugachev et al 2000, 2001, Pugachev and Xing 2001, 2002). In their earliest proposal, 
a non-iterative filtered back projection (FBP) method, using simple exponential 
attenuation for parallel X-ray beams in their dose calculation, was used to solve the 
inverse problem in IMRT while the beam orientations were optimized using a simulated 
annealing method. This fast FBP approach, however, seems to preclude the use of OFs 
other than a dose-based OF. In Pugachev et al (2001), an iterative method (Xing and 
Chen 1996) was used to optimize the beam profiles and non-coplanar beams were used 
in the search space. A typical prostate case in this study required about 100 hours to be 
optimized by searching over 5000 combinations of beam angles. The use of the beam's- 
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eye-view dosimetrics (BEVD) (Pugachev and Xing 2001, 2002) allowed them a more 
efficient sampling of the search space: The sets of beams that have low BEVD score, 
i.e. those that are not likely to produce an optimal dose distribution, are immediately 
rejected without going through the beam-profile optimization. It was claimed that the 
BEVD-guided sampling can improve the speed by a factor of 10 while maintaining the 
same capacity to locate the global minimum as that of an unguided sampling. 

It is interesting to note that most of previous approaches to BAO used the dose- 
based OF (Bortfeld and Schlegel 1993, Rowbottom et al 1998, 1999a, 1999b, Pugachev 
2000, 2001, Pugachev and Xing 2002). Das et al (2003) proposed an IMRT beam 
orientation selection method that is based on target equivalent uniform dose (EUD) 
maximization. The objective is to maximize the EUD on the target while simultaneously 
imposing DV constraints on the OARs. These two objectives were combined into a single 
OF by a simple scale factor and it would be interesting to compare the performance of 
this combined OF with that of a strictly DV-based or EUD-based OF. 

The approach that we take in this paper is to use a fast dose calculation engine 
that significantly reduces the time required for the IMRT optimization of each given 
set of beam angles. This algorithm relies on the use of a compact, albeit approximate, 
dose kernel (Wu et al 2003). The main purpose of this paper is to show that this 
fast algorithm is sufficiently accurate to be used as a dose engine in a BAO. Although 
simulated annealing has been used in several previous works on BAO, in this paper 
we show that the combination of simulated annealing, in particular the FSA algorithm, 
and the fast dose engine that we proposed recently can provide a practical, i.e. fast 
and useful, BAO system. In addition, by focusing the improvement on the speed of 
the dose engine in our system, we are also able to retain all other aspects of normal, 
i.e. clinical, IMRT optimization in our BAO. Thus our BAO system allows a DV-based 
OF, or indeed any other types of OF, to be used for the IMRT optimization for each 
given set of beam angles generated during the BAO. The promise of a BAO system is 
to provide better alternatives to plans produced by a planner, which in clinical settings 
are commonly generated by a DV-based optimization. It is therefore important to use 
a DV-based, instead of a dose-based, optimization in a clinical IMRT BAO system so 
that a direct comparison can be effected. If the BAO plan is deemed by the planner 
to be significantly better, it can be used for treatment; otherwise, the plan produced 
manually by the planner is retained. Thus the BAO can be used to automatically check 
if the best plan of the planner can still be improved by beam-angle adjustment. Starting 
with an already acceptable plan, the BAO is executed with the same DV constraints 
and parameters and its plan will be used if it is significantly better and can be produced 
within an allowable time. 

It should be noted that no a priori assumption is made regarding the selection 
of beam orientations in this work. Although incorporation of some criteria, such as 
constraints on the angular separation between adjacent beams, is likely to improve the 
performance of the system, we feel that such topics require separate careful studies. 
These additional beam-selection criteria will be incorporated in further studies. 
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Figure 1. Flow diagram of a BAO using the simulated annealing method. The system 
has been implemented to accept non-coplanar beams in the search pool shown in the 
topmost box. The kernel is calculated using an accurate SC dose engine of the TPS 
for each beam in the search pool. This kernel is used during the IMRT beam-profile 
optimization for dose calculation using the TLP method. Note that in the final step, 
once the optimal set of beams is selected, a final IMRT beam-profile optimization is 
performed using the SC method before the plan is output. 

2. Methods and materials 

In this work, we propose an efficient BAO method that combines a fast dose calculation 
engine, which we call the table lookup (TLP) method (Wu et al 2003), with an effective 
global search method based on the fast simulated annealing (FSA) algorithm. We have 
recently implemented this method in a program that interfaces with our in-house clinical 
IMRT system and studied its performance and potential. 

The following subsections will detail the components of our BAO implementation. 
Section l2?T] discusses the flow diagram of our BAO system and the algorithm that we use 
to select successive sets of beams from the search pool (the FSA algorithm). The dose 
engine that we use as a workhorse in the BAO (the TLP method) is presented in section 
12. 21 while section EHl describes the DV OF that we use to guide the IMRT optimization. 
Details on how the TLP method is used to accelerate our IMRT optimization and 
discussion on its accuracy are provided in section 12.41 The cases that were used to test 
our BAO system in this work are outlined in section 12.51 

2.1. Beam angle selection with fast simulated annealing 

Figure Q displays the flow diagram of our BAO system which is interfaced with an in- 
house IMRT system. The IMRT system, in which the planner sets general optimization 
parameters including the number and energy of beams whose directions are to be 
optimized, is in turn interfaced with a commercial treatment planning system (TPS) 
Pinnacle^ (Philips Oncology System, CA, USA). This in-house IMRT BAO program, 
which is external to the TPS, communicates with the TPS through an application 
programming interface, provided by the TPS vendor, called PinnComm. This allows 
the BAO and the IMRT programs to, e.g., request the TPS to perform a dose calculation. 

The BAO is performed using an FSA algorithm. Prior to entering the SA loop, the 
planner deflnes a search pool of beams. We use 180 equispaced beams for a coplanar 
search pool and 216 beams (72 equispaced beams at couch angle of 0°, 45°, and 315°) 
for a non-coplanar search in this work. For each beam in the search pool, the TPS is 
called by the BAO program to calculate an accurate dose for a uniform intensity profile. 
The jaws are automatically set for each beam to sufficiently cover the target with some 
margins. The kernel is then extracted from the TPS using the TLP method (section 
12. 2|) . Note that we do not optimize the number of beams in this study. 
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2.1.1. Fast simulated annealing: temperature lowering schedule. The BAO selects the 
optimal set of beams from the search pool. Successive sets of beam angles are generated 
during the BAO by using the FSA algorithm (Szu and Hartley 1987). The efficacy of 
this global search method has been proven in many fields, including IMRT and beam- 
angle optimization (Webb 1989, 1991, Mageras and Mohan 1993, Xiang and Gong 2000, 
Pugachev 2000). The FSA is known to be much more efficient than the classical SA 
(CSA) algorithm (Kirkpatrick et al 1983). To converge to a global minimum with 
probability one, the CSA requires the annealing temperature to be lowered with an 
inverse-logarithmic schedule: T = Tq/ ln(l + i), where i is the iteration number (Geman 
and Geman 1984). On the other hand, the FSA allows the temperature to be lowered 
with a much faster schedule (Szu and Hartley 1987): 



The SA temperature T is used in determining whether a new set of beam angles is to 
be accepted or rejected. Note that temperature in this context should be understood 
only as an optimization parameter and it does not have any relation, except in analogy, 
to the temperature of any physical system. A move generated by the FSA is accepted 
with a Metropolis probability (Metropolis et al 1953): 



where AE" is the change of the OF (score). Thus a downhill move is always accepted 
while an uphill move is accepted with a probability e~^^l'^ . The Metropolis probability 
allows the optimization to simulate the behavior of a physical system in a thermal 
equilibrium while the SA temperature schedule controls the lowering of temperature 
such that the system is effectively at a thermal equilibrium at any time. As in a 
physical system, this recipe ensures that local minima are properly sampled and the 
global minimum is reached as the temperature is lowered to zero. 

It should be noted that !\E in equation|2is the difference between the optimal score 
(final value of objective function) of the current trial set of beams and that of the last 
accepted set of beams. Each of these scores is the output of a beam-profile optimization 
performed on the corresponding set of beams. Details of the beam-profile optimization 
is not needed by the FSA algorithm whose task is to generate subsequent sets of beams 
and accept or reject each set based on the value of its optimal score. 

2.1.2. Fast simulated annealing: sampling distribution. The relative efficiency of the 
FSA over the CSA is related to the sampling distribution that it uses. At each BAO 
iteration, FSA samples the changes to the current beam angles x = (^i, ^2? ■ ■ ■ ? ^n) from 
a long-ranged Cauchy-Lorentz distribution: 



where n is the number of beam angles (the dimension of the vector x). In contrast 
to FSA, the CSA uses a short-ranged Gaussian sampling distribution: p(Ax) oc 



p = min [l,e-^^/^]. 



(2) 



p(Ax) oc [(Ax)' + W^(T)2]-("+^)/2 



(3) 



Algorithm of a clinical beam-angle optimization system 



7 



exp[— (Ax)^/2W^(T)^], where W{T) is the temperature-dependent width of the samphng 
distribution. Search-space regions farther away than W (T) from the current position 
are practically not sampled by the CSA distribution and therefore CSA requires a larger 
number of iterations than FSA to thoroughly sample the entire search space. 

Several different choices for the definition of W{T) are available in the literature 
(Mageras and Mohan 1993) and, in general, the "best" choice for an application can 
only be found by trial and error. For the FSA work reported in this paper, we have used 
a simple form W{Ty — oT"^. The value of the range parameter a was chosen such that 
even at the lowest temperature (T = 1 in this paper) each beam direction still has 1% 
probability to jump to another direction "farthest away" from it (e.g., from 0° to 180° 
gantry angle). Thus at any temperature, the current trial set of beams always has a 
finite probability to jump to any other point in the search space. The ability to sample 
regions beyond the local neighborhood is imperative in the simulated annealing method 
to avoid being trapped in a local minimum. 

In generating the random move for the set of beam angles x = {9i,92, - ■ ■ , 9n) we use 
the approach where the orientations of all beams in the current set are simultaneously 
updated. This has been shown (Mageras and Mohan 1993) to be more effective than 
the alternative approach, in which one beam is picked randomly from the current set 
and its orientation is moved by a random distance which is picked from the sampling 
distribution (Das et al 2003). 

2.2. Fast dose-calculation method 

The total dose at a given voxel, Di, is related linearly to the beamlet weights Wj by the 
dose kernel K^f 

D,^J2KijWj. (4) 
j 

Note that in this equation the beamlet index j comprises also the beam index; the index j 
runs over all beamlcts in all beams. The dose kernel Kij is a sparse matrix. At each dose 
voxel, significant contributions to its dose come only from a limited number of beamlets 
in each beam, the largest contribution being from the beamlet that directly traverses 
the dose voxel. The sparse nature of the kernel matrix can be shown numerically by, 
e.g., plotting the magnitude of the elements of the square matrix K^K (Crooks and 
Xing 2001). This also gives a measure of the overlap between beamlets. 

The fact that the kernel matrix is sparse has been used to devise a fast, albeit 
approximate, table-lookup (TLP) dose calculation method, that is useful for IMRT 
optimization (Wu et al 2003). This efficient method was proposed to improve the 
performance of applications that involve a great number of dose calculations. In the 
present work we have used this method as the fast dose calculation engine for our 
IMRT BAO. In the TLP method, each dose voxel is assigned one, and only one, 
beamlet from each beam, namely the beamlet whose raytrace intersects the voxel. 
The dose to the voxel is initially calculated using an accurate method, e.g., the 
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superposition/convolution (SC) method, for a uniform intensity profile for all beams. 
The TLP kernel is then calculated using the following formula: 



Once the kernel is calculated, dose calculation in subsequent iterations is done by looking 
up the kernel. Note that in addition to improving the speed of the calculation, the TLP 
method also greatly reduces the number of kernel elements that need to be saved and 
read, and therefore the time needed to access them. This number, for each beam, is equal 
to the number of dose voxels because of the one-to-one dose-beamlet correspondence (Wu 
et al 2003). The use of this fast dose-calculation method in our IMRT optimization will 
be further described in section 12.41 

2. 3. Objective function 

For each trial set of beam angles generated during the BAO, an IMRT optimization is 
performed to optimize the beam profile. This IMRT optimization is implemented by 
using an objective function, or score, that provides a measure of "distance" between 
the current dose distribution and the desired distribution. Starting from an initial 
distribution, the objective function is iteratively minimized by the optimizer. 

Several types of objective function (OF) exist in the literature and implemented 
in clinical IMRT optimization systems. Among these are the dose-based, DV-based, 
and the biology-based functions. The dose-based OF is the most straightforward 
conceptually and the easiest to implement (Webb 1989): it uses a penalty which is 
a quadratic function of the difference between the actual dose and the desired dose, 
which is usually the prescription dose for the tumor and zero for the organs at risk 



The DV-based objective function is among the ones that are most useful in IMRT 
optimization (Bortfeld et al 1997, Spirou and Chui 1998, Wu 2000). It has the advantage 
of being more intuitive than the biology-based OFs and more flexible than, since it is 
the superset of, the dose-based OF. The DV-based optimization typically requires more 
iterations to achieve convergence than the dose-based optimization. Furthermore, it 
requires the dose-volume histogram (DVH) to be computed after each dose calculation 
which adds to the total optimization time. In general, however, DV objective requires 
less adjustment of the optimization parameters to achieve a suitable plan (Wu 2000). 
Although multiple local minima exist in a DV-based IMRT optimization (they are absent 
in dose-based optimization), it has been shown (Wu and Mohan 2002, Llacer et al 
2003) that they are not expected to affect the outcome of optimization using gradient 
techniques in any clinically significant way. 

For the purpose of establishing our notation for later reference in the rest of this 
paper, here we present a brief description of our implementation of the DV OF. There 
are two types of constraint in a DV-based optimization, each is specified by a set of 




if voxel i is on the path of beamlet j, 
otherwise. 



(5) 



(OARs). 
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dose- volume values {Da, H)- The first type of constraint (upper constraint) requires the 
volume receiving dose greater than Da, which is denoted by V{Da) = Va, to be less 
than Vb- This constraint is implemented by searching for the dose value Df, such that 
V{Di,) = Vb- A penalty is charged if Va > Vb, which is equivalent to Da < Db since 
the DVH is a non-increasing function of dose. The upper constraint for a DV-based 
objective function therefore has the following form: 

/i = pi ^ H{Db - A) • H{Di - Da) ■ (A - Daf, upper constraint. (6) 

i 

Here H{x) is a Heaviside step function defined by 

fl, x>0. 




An upper constraint is used to limit the dose to an OAR or the hot-spot regions in a 
target. The index i in the summation runs over all voxels in the OAR, while pi is the 
penalty prefactor. The constraint's objective is to reduce the volume of the OAR that 
receives dose greater than Da, it therefore penalizes voxels that have dose Di > Da- 
Note, however, that voxels that have dose greater than Db are not penalized since the 
total volume of such voxels is already lower than Vb- 

A lower constraint is used to limit the cold-spot regions in a target. This constraint 
requires Va, the volume receiving dose greater than Da, to be greater than Vf,. The 
penalty is imposed if < Vf,, or equivalently if Da > Db'. 

/2 = P2 E ^(^a - A) • H{Di - Db) ■ (A - A)', lower constraint. (8) 

i 

In an IMRT optimization, DV constraints are assigned by the planner to the 
volumes of interest (VOIs). Each VOI can have more than one constraint; these 
constraints control the shape of the DVH curve for the corresponding VOI. Since there 
is no specific desired dose for each VOI in a DV-based optimization, in contrast to 
dose-based optimization, the quahty of the DVH of one or more VOIs can be sacrificed 
slightly by the optimizer to improve those of some other VOIs. The importance of each 
DV constraint, relative to the other existing constraints, is set by the relative value of 
its penalty prefactor. 



2.4- Intensity optimization for each set of beam angles 

For each trial set of beam angles in the BAO, an iterative IMRT optimization using the 
Newton's method is performed to optimize the beam profiles by minimizing a DV-based 
objective function. At each iteration during the IMRT optimization each bcamlct weight 
Wj is updated to Wj + SWj. The size of this update 6Wj depends on the current values 
of dose and the kernel: SWj — F[{Di}, {Kij}]. The procedure for computing 6Wj can be 
found in Wu and Mohan (2000) and will not be repeated here. Our IMRT optimization 
program is interfaced with a commercial 3D TPS. For a given set of beams, we start 
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from a uniform intensity profile for all beams and compute an accurate dose {-Dj} using 
the dose calculation engine (which is based on the SC method) of the TPS. Kernel {Kij} 
is then calculated for the iterative optimization using the TLP method (equation • 

2.4-1 ■ Intensity Optimization with Superposition Convolution (lO-SC) Method. Note 
that dose and kernel enter independently in the computation of the beamlet weight 
update 5Wj. This allows different levels of approximation to be used for dose and 
the kernel. Ideally both dose and kernel are computed accurately using, e.g., an SC 
method. However, for an IMRT system that is interfaced with a commercial TPS, the 
kernel is usually not directly available from the TPS, although the dose is. In the 
lO-SC method (Wu and Mohan 2000), 5Wj is obtained using an accurate dose and 
an approximate kernel at each iteration. The accurate dose, that incorporates scatter 
effects, is obtained from the SC dose engine of the TPS. The approximate kernel, on the 
other hand, is obtained from the accurate SC dose using equation El Although scatter 
effects are neglected in this approximate kernel, they are partially accounted for in the 
calculation of 5Wj since a full-scatter dose is used. The TLP method of constructing 
an approximate kernel has been found to be useful in practice (Wu et al 2003) and this 
lO-SC method has been used for IMRT routinely in our clinic. The main limitation 
of the lO-SC method is the fact that the dose has to be computed at each iteration 
using the SC dose engine from the TPS. The total optimization time is then strongly 
dominated by the lengthy time spent on repeated dose calculations. 

2.4.2. Intensity Optimization with Table Lookup (lO-TLP) Method. For the purposes 
of BAO, we introduce a further approximation to the computation of 5Wj which provides 
a significant speed up in the calculation time. Two different methods of dose calculation 
are used in the lO-TLP method. An initial dose is calculated accurately using the SC 
method of the TPS. A kernel is then generated from this accurate dose using equation 
El In subsequent iterations, instead of computing a fully accurate dose using the SC 
method, dose is calculated by looking up the kernel, i.e., using equation jH Thus both 
dose and kernel are approximate in the lO-TLP method as opposed to only the kernel 
being approximate in the lO-SC method. Note that scatter effects are not totally absent 
in the lO-TLP method because a full-scatter dose is always used as the initial dose for 
the beam-profile optimization. 

The TLP method provides a substantial speed advantage over the SC method. 
Since it is based mostly on memory lookup operation, each TLP dose calculation, which 
is the dose engine in the lO-TLP method, is typically at least two orders of magnitude 
faster than the corresponding SC calculation, which is used in the lO-SC method. 
In comparing the speed of these two optimization methods, we need to consider two 
separate situations: (1) a single IMRT optimization for a given set of beams; and (2) 
the BAO in which a large number of sets of beams (1000 in this work) are optimized 
consecutively. In the single IMRT case, the lO-TLP method uses the TLP dose engine 
to calculate the total dose at each iteration except the first one where SC is used to 
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Figure 2. Good correlation exists between tlie optimal IMRT score obtained using 
the lO-TLP method and that obtained using the lO-SC method (see main text for 
description of these methods). 

compute the kernel. The lO-SC method, on the other hand, uses the SC dose engine at 
each iteration. The speed ratio between the lO-TLP and lO-SC methods in this case 
is given roughly by the total number of iterations needed to converge the optimization, 
which is usually about 10-20 in our IMRT system. In the case of BAO, the kernel for 
each beam in the search pool is computed, using the SC method, prior to any IMRT 
optimization. During the beam-profile optimization of each given set of beams, no 
SC dose calculation is needed because the kernels have all been computed and can be 
quickly read from memory by the TLP method. Except for the time overhead spent in 
the initial calculation of the kernels, which is small relative to the total BAO time (table 
12} the speed ratio in this case is given directly by the speed ratio between the TLP and 
the SC methods. Thus our BAO system is at least two orders of magnitude faster than 
it would be if we used the SC method as the dose engine for beam-profile optimization. 

The feasibility of using this approximation for a BAO is demonstrated in Fig. 
12] where we compare the optimal scores obtained using the lO-SC method with the 
corresponding scores obtained using the lO-TLP method. This plot has been generated 
using 100 randomly selected sets of 3 beams that were used to plan a typical prostate 
case. Each score shown in Fig. |21 is the square root of the optimal value (according 
to the respective optimization method) of the objective function which is the sum of 
contributions from all DV constraints applied to the VOIs in this case (Wu and Mohan 
2000). It can be seen that the lO-TLP scores are well correlated with the lO-SC scores: 
"good" or "bad" sets of beam angles according to the lO-SC method are predicted 
similarly by the lO-TLP method. In general the lO-TLP scores are lower than the 
lO-SC scores. This is due to the fact that less scatter components are included in the 
lO-TLP method than in the lO-SC method. Dose "leakage" outside the target, due 
to scatter, is slightly underestimated by the lO-TLP method and this gives rise to the 
undercounting of the score and the finite variance of the correlation. It should be noted 
that the systematic score undercounting in the lO-TLP method is practically irrelevant 
to a BAO since only the relative ranking of different sets of beams is needed. The 
presence of a finite variance in the correlation indeed introduces some limitation to the 
predicting power of the lO-TLP method. It is the purpose of this paper, however, to 
show that even with this limitation the lO-TLP method can still be used as an effective, 
and efficient, IMRT workhorse in a BAO system. 

It can be seen in figure ^ that our BAO algorithm simply serves as an additional 
"outer loop" to a normal IMRT optimization (Wu and Mohan 2000) . The beam-profile 
optimization component of the BAO, however, uses a different dose-calculation engine 
than the one that is used in normal IMRT optimization. The lO-TLP method is used 
as the beam-profile optimizer during a BAO, because of its speed, while a normal IMRT 
optimization (for a given set of fixed beam angles) should be performed using the lO-SC 
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Figure 3. Typical progression of the accepted score in a BAO using an FSA algorithm 
as a function of the iteration number. The initial score is set to Tq — 1000 in this run. 
This is more than enough to allow all possible SA moves. As the optimization proceeds, 
the temperature is lowered with a To/(l + ?) schedule. As the temperature is lowered, 
SA practically turns into a local search algorithm although in FSA any point in the 
search space always has a finite probability to be sampled at any temperature. 

Figure 4. DVH comparison for coplanar beams, (a) Significant improvement can be 
obtained for a 3-beam prostate plan. The DVHs for 3 equispaced beams are shown 
by the dotted and dash-dotted lines, while those for the optimal BAO beams are 
shown with solid lines. BAO provides clear improvement for all of the DVHs. (b) 
Comparison of different sets of equispaced beams. In this case, improved sparing of 
the rectum results in worse sparing of the bladder, or vice versa. The DVHs have been 
normalized to 72 Gy at 90% of target. 



Table 1. Comparison of the prostate and the head-and-neck cases that have been 
used in this work to test our BAO system. 

Parameter Prostate Head-and-neck 



Number of CT slices 

Dose resolution 

Dose matrix dimension 

Size of bitmap vector 

Intensity matrix size (per beam) 

Intensity resolution 

Target volume (cc) 

OAR volume (cc) 



26 

0.2 X 0.2 X 0.2 cm^ 

(51,87,62) 

98,985 
- (20 X 15) 
0.3 X 0.5 cm^ 
91 (Prostate) 
252 (Bladder) 
78 (Rectum) 



129 

0.2 X 0.2 X 0.2 cm^ 
(175,91,119) 

388,238 
- (85 X 45) 
0.3 X 0.5 cm^ 
62 (GTV), 175 (CTV) 
24 (Cord), 843 (Nodes) 
32, 17 (Right, Left Parotid) 



method, because of its accuracy. 

It should be noted that although the lO-TLP method is used during the BAO, a 
final beam-profile optimization using the lO-SC method is performed on the optimal set 
of beam angles output by the BAO. This is a part of the finalization step shown in the 
last box in figure H This final optimization of the BAO plan facilitates its comparison 
with a human-prepared plan, which is also optimized using the lO-SC method. In the 
following section, all of the displayed results have been produced using a lO-SC-based 
optimization. 



2.5. Materials 

Two cases have been used in this work to test the practicality of our BAO system: a 
prostate case and a head-and-neck case. Some of the system parameters of these two 
cases are displayed in table [U The prostate case represents the cases with relatively 
small size and simple arrangement of the target and the OARs. The head-and-neck 
(HN) case represents the other end of the spectrum namely those with relatively large 
size and with multiple targets that are located close to sensitive OARs. For the prostate 
case we compare the BAO results for 3, 5, and 7 coplanar beams with that of the 
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Figure 5. DVHs produced by BAO for (a) 5; and (b) 7 coplanar beams compared 
to those produced by the corresponding number of equispaced beams for a prostate 
case (the same case as in figure 0)1 . The non-equispaced beams selected by the BAO 
are better able to satisfy the specified DV criteria (table OJ. The DVHs have been 
normalized to 72 Gy at 90% of target. 



corresponding number of equispaced coplanar beams. We also examine the dependence 
of the improvements that can be obtained from BAO on the criteria used to guide the 
optimization. This is done by studying the results of two different sets of criteria: a 
"soft" and a "hard" set. For the HN case we compare the BAO results for 7 and 9 
non-coplanar beams with that of 9 equispaced coplanar beams. 

The computational time required to perform a BAO on a given case generally scales 
with the system size. This is determined principally by the product of the total number 
of beamlet intensities that need to be optimized and the total number of voxels for which 
dose needs to be computed during the optimization. Note that table [T] lists the typical 
size of the intensity matrix (which is the size of the smallest rectangle that encloses the 
target in the BEV of the beam) since the specific size generally varies with the orientation 
of the beam. Furthermore, the number of beamlet intensities that are optimized is 
generally less than the size of the intensity matrix because the beamlets that do not 
pass through any voxel in the target can be completely turned off. Although increasing 
the number of optimization criteria is expected to increase the total optimization time, 
the increase is minimal because once the dose is calculated the value of the OFs can 
be computed very quickly. The size of the bitmap vector for the HN case in table ^ is 
about 4 times that of the prostate case and the size of its intensity matrix is about 4x3 
times larger. Thus, for the same number of beams, the total IMRT optimization time 
for the HN case is expected to be about 48 times that of the prostate case. 

3. Results 

3.1. Prostate case 

Figures and El display the improvements that can be obtained from BAO for a prostate 
case. Each of the results shown in these figures has been obtained from a 1000-iteration 
FSA run and was completed in about 1-2 hours on a 400 MHz workstation (table 
121). The initial SA temperature Tq was set to an arbitrary large number (Tq = 1000) 
which allowed practically all random moves generated by the SA sampling distribution 
to be accepted. This is an important practical requirement of the SA algorithm: The 
initial temperature has to be high enough to allow random sampling of the search 
space during the initial iterations. Subsequently, the temperature was lowered using a 
To/(l +i) schedule. Figure El displays a typical progression of the score with increasing 
iteration number. The optimization is stopped after 1000 iterations (final temperature 
is 1) because practically the temperature is already too low to allow a further significant 
change in the score to occur. Note that, based on our experience as shown in figure 
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Table 2. Comparison of user CPU times for BAO with 1000 sets of beams. Note 
that additional time is also needed to calculate the kernel for the beams in the search 
space: about 15 minutes for 180 beams in the prostate case and 2 hours for 180 for 



216 beams in the HN case. 


Case 


CPU Speed 


User CPU Time 


3-Beam Prostate 


400 MHz 


1 h 5 min 


5-Beam Prostate 


400 MHz 


1 h 40 min 


7-Beam Prostate 


400 MHz 


2 h 20 min 


7-Beam HN 


750 MHz 


13 h 15 min 


9-Beam HN 


750 MHz 


15 h min 



Table 3. DV criteria used for optimizing the prostate plan. 



VOI 


Type 


Dose (cGy) 


Volume (%) 


Penalty 


Target 


lower 


7200 


90 


10 




upper 


7900 


5 


3 


Bladder 


upper 


6000 


1 


5 




upper 


4000 


5 


3 


Rectum 


upper 


6000 


5 


5 




upper 


4000 


8 


5 



El the optimization can generally be terminated after a few hundred iterations without 
incurring a significant cost in the final score. 

The optimal beams for this prostate case have been chosen from 180 coplanar 18- 
MV photon beams (at zero couch angle) which are distributed uniformly over the 360° 
range for the gantry angle with a 2° interval. To facilitate the comparison, each plan 
shown has been normalized to 72 Gy at 90% of the target volume. The DV-criteria that 
were used for the IMRT optimization are listed in table El 

Figure compares the DVHs of a set of 3 equispaced beams and those of a set 
of non-equispaced beams selected by the BAO. It can be seen that the BAO beams 
provide significant improvements over the equispaced beams. The target receives a more 
uniform coverage while both of the OARs, the bladder and rectum, receive significantly 
less dose. Since in a 3-beam equispaced set each beam is separated from its neighbors 
by a relatively large 120° interval, it is important to examine if better results can be 
obtained from equispaced beams by uniformly rotating the complete set. This is done 
in figure Et) where we compare 4 different sets of equispaced beams, each equispaced 
configuration is rotated by 30° from the previous one. While their target coverages are 
similar, they spare the OARs rather differently. It is interesting to observe that the set 
that provides better protection for the bladder does so at the expense of the sparing 
of the rectum, or vice versa. In addition, all of these plans are still inferior to that of 
the BAO beams shown in figure E^. These comparisons illustrate the importance of 
non-equispaced beam- angle selection in IMRT. 
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Table 4. BAO gantry and couch angles (in degrees) for noncoplanar HN plans. 



Beam Number 


1 


2 


3 


4 


5 


6 


7 8 


9 


Gantry (9 Beams) 


5 


100 


140 


205 


240 


245 


290 300 


345 


Couch (9 Beams) 





45 


45 


45 


45 


315 


45 45 





Gantry (7 Beams) 


20 


105 


115 


245 


320 


325 


340 




Couch (7 Beams) 


315 





45 





45 


315 


45 





Table 5. DV criteria used for optimizing the head-and-neck plan. 



VOI 


Type 


Dose (cGy) 


Volume (%) 


Penalty 


GTV 


lower 


7300 


99 


50 




upper 


7700 


0.1 


30 


CTV 


lower 


6200 


96 


50 




upper 


6900 


40 


20 


Nodes 


lower 


5600 


95 


50 




upper 


5800 


15 


20 


Cord 


upper 


3000 


0.1 


50 


Brainstem 


upper 


3000 


0.1 


50 


Left Parotid 


upper 


2500 


30 


20 




upper 


5000 


1 


20 


Right Parotid 


upper 


1000 


32 


50 




upper 


3000 


10 


50 



Table 6. Prescription for the head-and-neck plan {R^ — 6810 cGy). 



DV Prescription (cGy) 



D98(GTV) 


> 6810 


D2(GTV) 


< 7491 


D95(CTV) 


> 6000 


Dgo (Nodes) 


> 5400 


Di(Cord) 


< 4500 


Di (Brainstem) 


< 5500 



DVH comparisons between equispaced and BAO plans for the same prostate case 
using 5 and 7 beams are presented in figures and respectively. These plans have 
been optimized using the same set of DV criteria as that used for the 3-beam plans in 
figure m Although it is expected that the importance of BAO will diminish as more 
beams are used (Stein et al 1997), these cases show that BAO can still provide rather 
significant improvements over the equispaced beams. We have also performed similar 
comparisons as in figure IHd between different sets of 5 and 7 equispaced beams. As is 
expected, the difference among different 5- and 7-beam sets is smaller than that for the 
3-beam set and it decreases as we increase the number of beams from 5 to 7. 
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Figure 6. Comparison between DVHs produced by 9 equispaced coplanar beams and 
those produced by 9 non-coplanar beams selected by BAO for a head-and-neck case. 
The DVHs have been normahzed to 60 Gy at 95% of the CTV volume. Solid lines are 
for the 9 non-coplanar beams selected by the BAO (tabled} and dashed lines are for 
the 9 equispaced coplanar beams. 

Figure 7. HN isodose contours for coplanar equispaced plan (top) compared with 
those of non-coplanar BAO plan (bottom). Nine beams are used in each plan. BAO 
plan provides significantly better sparing of the parotid glands compared to that of 
the equispaced beams. 

Figure 8. Comparison between DVHs produced by 9 equispaced coplanar beams and 
those produced by 7 non-coplanar beams selected by BAO for a hcad-and-neck case. 
The DVHs have been normalized to 60 Gy at 95% of the CTV volume. Solid lines are 
for the 7 non-coplanar beams selected by the BAO (tabled} and dashed lines are for 
the 9 equispaced coplanar beams. 

3.2. Head-and-neck case 

Due to its anatomical position, head-and-neck tumors are good candidates for the 
appHcation of non-coplanar IMRT. Limited experience with non-coplanar HN plans 
and the fact that selection of non-coplanar beam directions for IMRT is a highly non- 
intuitive process also make global search methods like the FSA, which is implemented 
in our BAO system, particularly attractive. 

Figures El and |H1 compare the DVHs for a HN case of 9 and 7 non-coplanar beams, 
respectively, with those of a benchmark 9 equispaced coplanar beams. The DV criteria 
used to optimize these HN plans are listed in table while the gantry and couch angles 
selected by our BAO algorithm for these sets are tabulated in table |^ Each of the 
BAO results for this HN case has been obtained from a 1000-iteration FSA runs. The 
optimal beams have been chosen from 216 possible directions: the search pool for this 
case consists of 3 possible couch angles (0°, 45°, and 315°) with 72 possible gantry 
angles distributed uniformly with 5° interval for each couch angle. In principle, finer 
resolution for the couch angle may be used. However, this will also increase the memory 
requirement since our BAO system, as it is currently implemented, requires the kernel for 
each beam to be available for quick access in the system memory. For this HN case, the 
size of the kernel is about 3 MB for each beam and therefore it requires about 650 MB of 
memory to store all the kernels. We have used a 750-MHz Sun workstation to perform 
the BAO for this HN case. The 1000-iteration BAO run for the 9 non-coplanar beams 
required 17 hours to complete (this includes about 2 hours of computing time required 
for kernel extraction) while the 1000-iteration run for the 7 non-coplanar beams required 
13 hours (kernels were read from disk for this case). The relatively long optimization 
time for this HN case, compared to the corresponding time for the prostate case, is due 
to the substantially larger size of the system. A side-by-side comparison of these two 
systems is presented in table [H Note that, as pointed out previously in the discussion 
regarding figure the number of iterations used can be substantially reduced without 
significant reduction in the optimal score. This has the potential of shortening the 
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Table 7. "Soft" DV criteria used for optimizing the prostate plans shown in figure 
The criteria for the target are the same as in table O while those for the bladder are 
"softer". Only one of the criteria for the rectum in tabled is used here. 



VOI Type Dose (cGy) Volume (%) Penalty 



Target 


upper 


7200 


90 


10 




lower 


7900 


5 


3 


Bladder 


lower 


6000 


3 


5 




lower 


4000 


10 


1 


Rectum 


lower 


6000 


5 


5 



Figure 9. Two sets of 3 beam angles that were selected by the BAO as "good" 
sets when "soft" DV criteria (table [T)) were used for the IMRT optimization. These 
illustrate the importance of good selection of DV criteria in BAO. The DVHs have 
been normalized to 72 Gy at 90% of the target. 

total optimization time by a factor of two. Further reduction in time can be obtained 
from using a more selective search, e.g. not allowing the beams to be too close to each 
other, inclusion of prior knowledge, and utilization of parallel processors. These possible 
avenues for improvements are currently being investigated. 

Figure IHl compares the DVHs for a set of 9 non-coplanar BAO beams with those of 
a 9-beam equispaced coplanar set. The BAO set provides both better dose uniformity 
at the GTV and CTV and better protection of the parotid glands. This better sparing 
of the parotid glands can also be seen in the isodose contours shown in figure [7| Figure 
IH] shows that a plan of similar quality can also be obtained using 7 non-coplanar beams, 
although in this case the quality of target coverage is lower than that of the 9 non- 
coplanar beams. The dose prescriptions for these plans are shown in table IHl The plans 
have been normalized to best satisfy all of these prescriptions. These have been achieved 
by normalizing the DVHs to 60 Gy at 95% of the CTV volume. 

4. Discussion 

4.I. Dependence on criteria 

A noteworthy feature of a DV-based optimization is that further improvement will no 
longer be attempted by the optimizer once the criteria are satisfied. In other words, 
once the contribution of a DV criterion to the OF reaches zero, it becomes irrelevant to 
the optimizer. Consider, e.g., a single criterion which is assigned to the rectum: < 6000 
cGy at 5% of its volume. Once this criterion is satisfied, e.g., when D5 (rectum) has 
been lowered to 5900 cGy, the optimizer will no longer attempts to improve the dose to 
the rectum, even though it is still possible to reduce its dose further. 

We have found that the improvement on plan quality that is output by our BAO 
system generally depends on the DV criteria that we assign to the VOIs. Figure IHl 
compares two of the top 5 plans generated by the BAO (for the prostate case presented 
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in section IXTj) using several "soft" criteria listed in table [7[ with the corresponding plan 
for a set of 3 equispaced beams. These plans have been generated for the same prostate 
case as that used in figures lU and The BAO sets are both predicted, by the lO-TLP 
method which is used during the BAO, to have lower score compared to the equispaced 
beams (55 and 75 for plans shown in figure IHK and IHb, respectively, compared to 186 
for the equispaced beams). However, upon reoptimization using the lO-SC method, the 
BAO sets produce higher score than the equispaced beams (186 and 212 for figure int- 
andEb, respectively, compared to 162 for the equispaced beams). 

The large score increase encountered in going from the lO-TLP to the lO-SC 
method for the plan shown in figure IHt is largely due to the additional score accumulated 
by the DVH of the rectum. In this case, the lO-TLP method minimized the score by 
creating a steep gradient of the DVH curve just below the 6000 cGy dose. Note that 
even though below 6000 cGy the BAO DVH for the rectum is higher, and therefore is 
worse, than the corresponding DVH of the equispaced beams, the contribution of this 
DVH to the OF can be smaller if its DVH at 6000 cGy is lower than the corresponding 
value for the equispaced beams. This is specifically what happens in this case due 
to the steep gradient of the BAO DVH curve. In going from the lO-TLP to the 10- 
SC method, however, this steep gradient cannot be sustained because of the inclusion 
of more scatter components in the lO-SC method. This "dose leakage" substantially 
increases the partial volume of the rectum that receives dose greater than 6000 cGy and, 
consequently, produces a large increase in score. 

The BAO plan in figure IHb satisfies the criteria listed in table [3 as well as, or indeed 
better than, that of the equispaced beams. In this case, the higher score of the BAO 
plan is due to the contribution from the "normal tissue" (NT) region. The NT region 
encompasses generic normal tissue surrounding the target that does not belong to any 
particular OAR (Wu and Mohan 2000). For the prostate case that we use in this work, 
the NT is a 1 cm extension of the target that is not part of any OAR and it is assigned 
an upper constraint of 6000 cGy at 5% of its volume. The "dose leakage" from the 
target to its surrounding NT, for the case shown in figure Eb is substantially larger than 
the allowance set by the assigned DV criterion and it therefore introduces a significant 
contribution to the OF on going from the lO-TLP to the lO-SC method. This is due 
to the small angular separation of two of the 3 beams in the BAO set (238° and 254°) 
that provides an elongated dose distribution along the direction of these beams. 

In contrast to this relatively large uncertainty in the prediction of the lO-TLP 
method for cases with "soft" DV criteria, in our experience we found that its prediction 
is much more reliable for cases with more demanding criteria, such as listed in table El 
and displayed in figures El and Although we only present this dependence on criteria 
using a prostate-case example in this paper, similar behavior has also been observed for 
the HN case. One possible reason that may be given for this observation is that the 
degeneracy of the solution, i.e. the multiminima of the OF, is larger for "soft" criteria 
(Wu and Mohan 2002). A careful study, however, is necessary in order to provide 
a definitive answer to this question. From the point of view of clinical applications. 
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the fact that BAO is more useful for "hard" cases is not expected to ehminate its 
practical relevance. It is well known that "soft" cases can be planned satisfactorily 
with equispaced beams or, indeed, class solutions. On the other hand, for "hard" cases, 
e.g. those of re-irradiation or of complicated geometry, class solutions are non-existent 
(almost by definition) and equispaced beams may not provide the optimal solution. For 
these cases the fast BAO described in this paper provides an attractive alternative to a 
manual trial-and-error approach. 

4.2. Further speed up? 

Several schemes have been proposed recently to rank beam directions according to their 
individual potential for creating a desirable IMRT distribution such as the single-beam 
cost function (Oldham et al 1998), beam's-eye-view dosimetrics (BEVD) (Pugachev 
and Xing 2001), and the figure-of-merit parameter (Chuang et al 2003). These schemes 
have also been studied as a possible way to intelligently accelerate a BAO. If a particular 
beam direction can be reliably predicted to always give rise to high values of objective 
function, independent of the directions of the other beams in the set it is a member of, 
it can be quickly removed from the search space. It is perhaps important to point out 
here, in connection with the OF of the optimization, that these schemes are expected 
to have less significance in a DV-based optimization than in a dose-based optimization 
for which they have been proposed and studied. Individual ranking of beam directions 
is useful in a dose-based optimization where each VOI is assigned a specific prescribed 
dose. In contrast, in a DV-based optimization each VOI is effectively assigned a specific 
range of acceptable doses. In this case, several different distributions of the total dose 
may have equally good value of the DV OF (in other words, the OF possesses multiple 
local minima and degenerate solutions). It is therefore less meaningful to attach a 
specific geometric "shape" , and therefore potentially "good" or "bad" beam directions 
for creating this shape, to the optimal distribution. 

No a priori knowledge has been used to guide the selection of beam angles in this 
work. The use of reliable prior knowledge is expected to further improve the speed and 
potential of our BAO system. However, the complexity of the BAO problem requires 
judicious selection, and therefore careful systematic study, of what "knowledge" is best 
to be used for each given case. As has been discussed in this paper, this systematic 
study is especially more pertinent to the DV-based IMRT optimization that we use in 
our BAO system, because of the presence of the multiple local minima of the OF, than 
to the dose-based optimization that is commonly used in other BAO systems. In our 
BAO system, the beam angles are selected from the preset values in the search space. 
Some prior knowledge can be incorporated in the initial selection of beams that are to 
be included in the search space, e.g., more and less beams can be positioned in favorable 
and unfavorable directions, respectively. Alternatively, the beam angles selected during 
the optimization can be restricted to have a minimum separation from each other. Prior 
knowledge can also be used to reduce the number of beams in the search pool. Although 
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large search pools are used in this study, 180 and 216 beams for prostate and HN case 
respectively, in a clinical setting one may not need to use such large numbers. Using 
a reasonable spacing of 10° between adjacent beams in a coplanar pool, e.g., allows us 
to reduce the size of the search pool to 36. This will significantly decrease the number 
of possible combinations of beam angles to be sampled during the BAO. In the current 
version of our BAO system, each beam in the trial set is given a uniform initial profile 
prior to the beam-profile optimization. Instead of this uniform intensity distribution, 
the initial beam profile can be set to a non-uniform one, chosen based on some prior 
knowledge. This has the potential of reducing the number of iterations needed for the 
beam-profile optimization and, consequently, the total BAO time. 

5. Conclusion 

In this paper we have demonstrated that a combination of a fast table-lookup dose 
calculation method and the FSA algorithm can be used in a practical BAO system. This 
method has the advantage of being able to produce clinically significant improvement 
within a reasonable time. We showed, e.g., that the beam angles of a 7-beam prostate 
plan can be optimized within about 2 hours on a moderate 400 MHz workstation. A 
DV-based OF is used to guide the IMRT optimization in our BAO system. This provides 
a direct link to our clinical IMRT system and facilitates the comparison between the 
BAO plans and those produced by human planners. Since the beam-angle optimizer in 
our system simply wraps around the IMRT beam-profile optimizer, other OFs that are 
available to our IMRT optimizer, such as EUD, can also be used for BAO in our system. 
Our BAO system can also search for the optimal set over non-coplanar beams. This 
feature may be particularly helpful to planners because of the insufficient experience 
with non-coplanar beams in treatment planning. 
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